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We consider a solvable model of the decay of scalar variance in a single-scale random velocity field. 
We show that if there is a separation between the flow scale fcg ow and the box size fccL, the decay 
rate A oc (febox/fcflow) 2 is determined by the turbulent diffusion of the box-scale mode. Exponential 
decay at the rate A is preceded by a transient powerlike decay (the total scalar variance ~ t~ 5 ' 2 if 
the Corrsin invariant is zero, t~ 3 / 2 otherwise) that lasts a time t ~ 1/A. Spectra are sharply peaked 
at k = fcbox- The box-scale peak acts as a slowly decaying source to a secondary peak at the flow 
scale. The variance spectrum at scales intermediate between the two peaks (feb ox k <C fcfl ow ) is 
~ k + ak 2 + . . . (a > 0). The mixing of the flow-scale modes by the random flow produces, for 
the case of large Peclet number, a k~ 1+s spectrum at k S> fcfiow, where 5 oc A is a small correction. 
Our solution thus elucidates the spectral make up of the "strange mode," combining small-scale 
structure and a decay law set by the largest scales. 

PACS numbers: 47.27.Qb, 47.10.+g, 05.40.-a, 95.10.Fh 
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I. INTRODUCTION 

The problem of the decay of passive-scalar variance 
has recently been reexamined in the literature follow- 
ing the realization that the decay rates, spectra, and 
higher-order statistics based on small-scale Lagrangian- 
stretching theories GI EJ IE are not consistent with 
either numerical [E IE IE IE] or experimental [E E3 re- 
sults in the long-time limit. Instead, the scalar decay 
is dominated by an eigenmodelike solution dubbed "the 
strange mode" p| because it combines intricate small- 
scale structure with globally determined decay rate and 
self-similar statistics (self-similarity is also seen in nu- 
merical simulations of the related problem of kinematic 
dynamo There has been a growing understand- 

ing 0, Il2l Il3l IbH Hll that the overall decay rate is set 
by the slowest-decaying system-scale modes. This brings 
to mind homogenization theory |l6j | , which considers the 
turbulent diffusion of passive scalar at scales much larger 
than the flow scale and where it is the largest-scale mode 
that decays most slowly. In this paper, we use a simple 
solvable example to demonstrate that the strange-mode 
decay rate is the rate of turbulent diffusion of the box- 
scale mode and show how the spectra of scalar variance 
accommodate both this box-scale diffusion and small- 
scale structure. 

Qualitatively, the key idea quantified by our theory is 
as follows. A scalar field whose variance is at the scale 
smaller than or equal to the scale of the ambient random 
flow is mixed at a rate determined by the Lyapunov ex- 
ponent of the flow — this is the Lagrangian-stretching 
approach. However, if the size of the box is larger than 
the scale of the flow, the scalar field can have variance at 



the scale of the box. The rate of transfer of this large- 
scale variance to the flow scale (turbulent diffusion) can 
be much smaller than the Lagrangian mixing rate, in 
which case this slow transfer sets the global decay rate. 

Our model emphasizes scale separation between the 
box and the flow. Our results are complementary to |17| , 
where the decay of a scalar field is studied with more 
generality (in two dimensions). 

We consider the advection-diffusion equation 



ft. 



u • V0 = t]A9, 



(1) 



with a random Gaussian white-in-time velocity field 
(u l {t,x)u j (t' ',x')) = S(t - t')K ij {x - x') known as the 
Kraichnan model [l8|. The mean scalar concentration 
has been subtracted — i.e., (6) = 0. For the Kraich- 
nan velocity, the angle-integrated scalar-variance spec- 
trum in d dimensions T(k) — J dil^ fc d ~ 1 (|£>(k)| 2 ) satis- 
fies an integro-differential equation valid at all k j34j: 



d t T{t,k) + (2rj + Ko)k 2 T(t, k) 



hkj J d d k'K i: >{k-k r )T(t, k'), 



(2) 



where K y (k) = k(/c)(<5 u — kikj/k 2 ) is the Fourier trans- 
form of K y (x — x') and kq = (1/ c?)k m (0) is the turbulent 
diffusivity (d is the dimension of space). 

In Sec.m we review the theory of scalar decay at small 
scales, which leads to the standard Lagrangian-stretching 
results. In Sec. IIIII the theory for a finite-scale flow is de- 
veloped — this is the main part of the paper. Concluding 
remarks are in Sec. IIVI 



II. SMALL-SCALE THEORY 
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If the Peclet number Pe = uk &ov/ /r] is large, 9 varies at 
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scales as small as Pe ' times the scale of the flow. This 
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small-scale structure can be considered in the approxima- 
tion of spatially linear velocity field [l^,l20l |. — viz., 



1 



/^'(y)~/^- ^2 [y 2 S ij 



1 



2/V 



(3) 



For the Kraichnan velocity, the Lagrangian-stretching 
theories Q, 0, 0, H amount to the approximation 10 • 
The scalar-variance spectrum satisfies a Fokker-Planck- 
type equation [l8Ll20| 



d * T=D iL 



2r]k 2 T, 



(4) 



where D = [(d — l)/2(d + 1)]k2- This equation is valid 
for fc 3> fcfiow in this limit, it either can be obtained 
from Eq. (J2J) by expanding the mode-coupling term on 
the right-hand side or derived directly by assuming linear 
velocity field [H. 

The solution of Eq. that decays at k — ► oo is 

T(t, k) = C x e-^k-^K^^ (k/k v ), (5) 

where C\ is a constant, K v (z) is the modified Bessel func- 
tion of the second kind, 70 = (d 2 /4)D, k n = (£>/2?7) 1/2 , 
and A is the nondimensionalized decay rate, which must 
be calculated by applying the correct boundary condi- 
tion at small k. If we assumed that the decay rate is 
fully determined by the small scales, a reasonable pro- 
cedure would be to choose some infrared cutoff fc* and 
require the flux of scalar variance through fc* [the square 
brackets in Eq. |0J] to vanish (cf. 21] ). This can be satis- 
fied only for A > 1. For fc* <C k v , the zero-flux condition 
becomes sm[(d/2)y/X — 1 ln(fc*/2fc, ) )] = 0. Placing the 
cutoff fc* at the largest zero, ensures that T(fc) is every- 
where positive. We get (cf. |ljj) 



A = 1 



{2-K/df 



[ln(fc*/2fc r) ) 



l + 0(l/log 2 Pe). (6) 



This implies a fc( d ~ 2 )/ 2 spectrum at k k v . If the scalar 
variance is initially at k 3> fcflow, these results (or their 
analogs for other model flows) hold during the initial 
stage of the scalar decay. However, in the long-time limit 
for cases in which the system (box) size is (several times) 
larger than the flow scale, both numerical simulations [l4[ 
and experimental results ^(j obtain much smaller decay 
rates and spectra with negative exponents. The conclu- 
sion is that the zero-flux boundary condition is incorrect 
and the decay rate A < 1 must be determined by match- 
ing the solution Q to the solution at nonlarge k where 
Eq. (@J is invalid and Eq. J2J) must be used instead. The 
spectrum at fcfl ow <C fc fc,j is then ~ k s ^ x \ where [from 
Eq. @] 



«(A) = -l + (rf/2)(l- v / T^A). 



(7) 



Note that for A -C 1, s(A) ~ — l + (d/4)A, which coincides 
with the formula proposed in . 



For the initial spectrum concentrated at k ~ k n , the 
period of validity of Eq. © is the time it takes the spec- 
trum to spread to k ~ fcflow The spreading can be 
shown to be exponentially fast with a rate ~ 70. Since 

1/2 

k v ~ Pe ' fcflow, the Lagrangian-stretching results are 
valid for t <C 7q logPe. 



III. THEORY FOR A FINITE-SCALE FLOW 

The challenge now is to find A by solving Eq. J5J. Let 
us specialize to three dimensions (theory in 2D is analo- 
gous) and choose a simple form for the velocity correlator: 
k(k) — N6(k — fcflow), where N — 15K2/167rfc| ow [note 
that K2 = (2/5)fcfl ow Ko]- This describes a Kraichnan en- 
semble of randomly oriented "eddies" of size fcflow- We 
set fcflow = 1 and carry out angle integrations in Eq. 
to get 

d t T(t,k) + {2 V + KQ )k 2 T{t,k) 



15 1, 

Ko k 
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|i-fc| 



~ k dk' 

— K(k,k')T(t 7 k'), (8) 



where K(k,k') = -fc' 4 + 2(l + fc 2 )fc' 2 - (1-fc 2 ) 2 . It is not 
hard to ascertain that Eq. JSJ reduces to Eq. when 
k 3> 1. Let us now consider the opposite limit fc <C 1 - 
i.e., the evolution of scale variance at scales much larger 
than the scale of the flow. In this limit, the integral in 
Eq. |(SJ| is dominated by the modes in the neighborhood 
of the flow scale fc = 1. Neglecting 77, we get 

d t T + k 2 T=^-k [ dq (fc 2 - q 2 )T(t, l + q) = S(t, fc),(9) 
4 J-k 

where time is rescaled iKofcf| ow — * t. The solution is 



T(t,k) = 



T(0,fc)+ / dt'S(t',k)e 



k A t' 



(10) 



In order to complete the solution, we must determine 
T(t, l + q). For q -C 1, it satisfies 



3 f 2+q dk' / ., fc' 4 
Q tT + T=- Ik' 2 - — 

4 J\ q \ fc' V 4 



T(t,k'). (11) 



As we shall see, the solution <|10[) is sharply peaked at 
fc = fcpeak ~ 1/v^- In an infinite system, this peak would 
move indefinitely towards ever smaller fc. In a finite 
system, fc pca k eventually becomes comparable to the in- 
verse system size. Strictly speaking, this means that one 
must solve the problem with a discrete set of modes and 
application-specific boundary conditions (cf. 0,E1II2])- 
Instead, we model the finite box by introducing an in- 
frared cutoff fcbox into our continuous theory. All lower 
integration limits in fc space are subject to this cutoff. 
This is not a rigorous operation, but it is a reasonable 
modeling choice as long as fc DO x <C 1. The time at which 
fcpcak ~ fcbox is t ~ Therefore, there will be 



3 



two asymptotic regimes of scalar decay: the transient 
stage t <C 1/^box' when unmodified continuous theory 
can be used, and the long-time limit t 3> l/fc 2 ox , when 
the box cutoff is important. Note that even for the tran- 
sient stage, we assume t 3> f; i.e., we consider times 
much longer than the "turnover time" of the flow. For 
spectra initially at small scales, we harden this condition 
to t log Pe, so that the Lagrangian-stretching theory 
ceases to be valid. 

Consider Eq. ljTU|l . Since fc <§C f , we can Taylor-expand 
the initial spectrum: T(0, fc) = C 2 k 2 + C 4 fc 4 + . . . . Here 
C 2 is the Corrsin invariant, C 2 cx J d 3 x (0(x)0(O)) [23|. 
Some aspects of the scalar decay differ for cases with 
C 2 > and C 2 = 0. We shall first develop our theory for 
C" 2 > 0. The case of C 2 = will be treated in Sec. lrTTDl 



Case of C 2 > 



Consider first the transient stage f -C t <C 
us assume that the dominant term in the solution (|f 0[1 is 



T(t, k) = C 2 k 2 e 



2-k z t 



(12) 



which peaks at fc pca k = 1/yt- We shall justify this as- 
sumption a posteriori. Let us now determine the flow- 
scale solution [Eq. (fill) ] . Assuming that the main contri- 
bution to the integral in Eq. (|f f[) is from t'<l (also to 
be verified later) and substituting Eq. I|f 2fl for T(t,k'), 
we get, to leading order in in 1/t (and neglecting 9jT), 



true that the time integral in Eq. (|f Ofl is ~ fc 4 , so the 
above estimates remain valid. 

Now consider -C k -C 1. In this limit, S{t, k) ~ 

(3/4)V7i : Ti(t)^ 1 / 2 fc 3 . Substituting into Eq. JJJJ), we get 



T(t,fc) 



C 2 fc 2 e~ fe * 



|V5F2i(*)*~ 1/2 A- (15) 



The second term becomes comparable to the first at 
k — (lni/i) 1 / 2 . At larger k (but still k < I), it replaces 
the solution i|12|) as the dominant asymptotic. The solu- 
tion HI5|) is uniformly valid across the transition region. 
Together with Eq. (|lip. this implies that the solution 113|) 
is valid for |g| < (lni/*) 1/2 . 

In the long-time limit t > l/fc^, Eq. (JTJJ is still 
valid. Again, we assume that the dominant solution is 
Eq. l(T2*|) . Its peak is at fc pca k = fcbox- In Eq. ljTT|) . the 
lower integration limit is adjusted to max{|q|, fcbox}, as 
explained above. The flow-scale peak is now confined to 
\q\ < fcbox- At these q, Eq. lfT3|) is replaced by 



T(t,l 



\ 3 C 2 , 2 



? 2 )e 



-feL 



(16) 



This solution gives the dominant contribution to S(t, x) 
[Eq. ©], so we have (integrating from —fcbox to fcbox) 
S(t,k) = (3/8)0^ k(k 2 - fc 2 ox /5)fcg ox exp(-fc 2 ox t), 
which we substitute into Eq. ltTU|l . The box mode obeys 



T(t, fcbox) 



^2 fcbox 



C2 fcbox 



lni 



(17) 



T(t, I + <?) 



3^ 

8 i 2 



-« 2t sri(«) f 



(13) 



where Ti(t) = T(t, I). This describes the neighborhood 
of the flow scale, where the coupling to the small-fc modes 
produces a secondary peak with the width ~ 1/ \ft. We 
shall see below that Eq. (|I3|I is, in fact, valid beyond the 
width of the peak and up to \q\ ~ (lni/i) 1 / 2 . We now 
substitute Eq. (JT^J back into Eq. JTDJl to see that Eq. Q 
is, indeed, the dominant solution. 

For fc < 1/v^, we get 5(t, fc) ~ fc 4 Ti(i). Equation iJTDJ 
becomes, to two leading orders in fc, 



T(t,k) 



C 2 k 2 



dt'T^t') k 



-krt 



(14) 



Since T\{£) decays faster than 1/t [Eq. I|13l) ]. its time in- 
tegral in Eq. (|14|l tends to a constant when t ^> 1 and is 
dominated by the initial stage of the evolution of Ti(i') 
[there is no divergence at t' = because the solution 113|) 
is only valid for t' 3> 1]. This time integral represents a fi- 
nite amount of scalar variance that is initially transferred 
from the flow scale to the large scales (small fc) . That the 
effect of nonlinear coupling only appears in the fc 4 term 
is a reflection of the conservation of the Corrsin invari- 
ant: the coefficient in front of fc 2 cannot be changed. We 
see that, as long as C 2 ^S> const x fc 2 ox , the first term in 
Eq. I|14|) dominates. Note that for fc ~ 1/Vt, it is still 



The time integral now has a logarithmic divergence, rep- 
resenting a small amount of transfer of scalar variance 
from the flow-scale mode to the box mode. This contribu- 
tion is not significant because the second term in Eq. I|17|) 
only exceeds the first at t = exp(10/3fc 4 ox ), which is un- 
physically large even for moderately small values of fcbox- 
The width of the box-scale peak, for which the decay 
law (|17|l is valid, is estimated by fc 2 — fc 2 ox < 1/t — i.e., 



< 



l/2fc box^- 



Outside the peak (fc 2 — k^ oy _)t » 1, we have 



T(t, fc) 



c 2 fcV 



k-l 



3 C 2 fc 2 - 
8 t fc 2 



fcoox/5 



fc 2 

K box 



fc a p- 

"'box 1 ' 



(18) 



The first and second terms are of the same order when 
fc 2 — fc 2 ox ~ \nt/t — i.e., fc — fcbox ~ lnt/fcbox^- In the 
intermediate scale range fcbox <C fc <C 1, we have 



T(f , fc) 



3 (h 
8 t 



^box e 



Ti(t)fc, 



box'' 



(19) 



Let us summarize what we have learned so far. We 
have been concerned with two narrow bands of modes: 
the flow-scale modes T(t, 1 + q), q <C 1, and the large- 
scale peak T(t, fc pea k) at fc pea k = l/v^j which became the 
box mode T(t, fcbox) in the long-time limit t 3> 1/fcbox- 
The width of the peak was ~ l/\fi when t <C l/fcb ox an( l 
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FIG. 1: (a) Scalar-variance decay: bold lines depict the total variance E(t), thin lines the flow-scale mode T(t, 1), and dashed 
lines the box mode T(t, fcbox). Dotted lines are theoretical slopes. Time is in units of % 1 . The C 2 = runs are the same as in 
Fig- HI (b) Effective decay rate X(t) = —d\xiE{t)/dt. in units of the Lagrangian-stretching decay rate 70 for various values of 
fcbox- Dotted lines correspond to A = (40/9)fc DOX . The initial spectrum was the same as for Fig. |2] 



~ l/^box^ when t ^> l/^box- The flow-scale modes could 
be assumed to be coupled solely to the peak because of 
the peak's sharp dominance of all other modes: indeed, 
T(t, fc pcak ) ~ tTx(t) > Tx(t). The width of the secondary 
peak at the flow scale was determined by fcpcak- This 
width can be parametrized by A{t) = (fc^ ox + 1/i) 1 / 2 . 

The large-scale peak and the flow-scale band are sin- 
gularities of the scalar-variance spectrum. Note that in 
the long-time limit t l/fc^ox' the- width of the flow- 
scale peak is A = fcbox, while the width of the box mode 
is 1/fcbox^ fcbox- In a finite system, the spacing of the 
modes cannot be smaller than fcbox, so it is, of course, 
unphysical to talk about variation of the spectrum at 
distances less than fcbox- In our continuous theory, the 
collapse of the singularities to profiles narrower than fcbox 
means that they should be interpreted as single modes. 
Formally, they can be represented as S functions: thus, 
for the flow-scale band 1 — A(t) < fc < 1 + A(t), we write 
T(t, k) = NiTi(t)A(t)5(k - 1), where iVi is a constant 
arising from integrating over the specific shape of the sin- 
gularity: A^i = ^pii in the transient stage and iVi =4/3 
in the long-time limit. 

B. Nonsingular spectrum 

Let us now consider the nonsingular part of the spec- 
trum. For all fc > A(t), we write 

T(t, fc) = NiA(t)Ti(t) [S(k - 1) + /(fc)] . (20) 



Substituting this decomposition into Eq. © and neglect- 
ing d\n[A(t)Ti(t)]/dt, we find the equation for /(fc): 

1 r 1+k dV 

where H(2~k) is the Heaviside step function [the term it 
multiplies comes from integrating 5{k — 1) and vanishes 
for fc > 2]. For fc <§; 1, we expand /(fc') under the integral 
around fc' = 1 and find that the lowest-order term is 
/(l)fc 2 . Clearly, /(l) > 0. Thus, the spectrum at scales 
intermediate between the box scale and the flow scale is 

T(t, fc) = mA^T^t) [(3/4)* + /(l)fc 2 + ...]. (22) 

Coefficients in higher-order terms will involve derivatives 
of / at fc = 1 . The first term in the expansion (1221) comes 
from coupling to the flow-scale mode and is consistent 
with the asymptotics (|15|l and l|19() . so the nonsingular 
solution connects smoothly to the large-scale peak. The 
interaction of the nonsingular modes between themselves 
enters in the second order. 

To complete the solution, we would have to find /(fc) 
at fc > 1. For fc < 2, this means solving the full in- 
tegral equation, but we do not really need to do this. 
Note that the modes with fc > 2 are not directly cou- 
pled to the flow-scale mode. In a rough way, it can be 
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said that the Fokker-Planck regime starts at k = 2 and 
the solution of Eq. @ must be matched to T(t, 2) = 
N 1 A(t)T 1 (t)f(2). The specific value of /(2) is not impor- 
tant. The matching is done by setting C\ cx K{t)Ti{t) 
and A = — 7 ~ 1 d\n[A(t)Ti(t)]/dt (the effective decay rate) 
in Eq. Q. Here 70 = (9/16)« 2 = (9/40)fcg ow K is the 
decay rate from the Lagrangian-stretching theory. Dur- 
ing the transient stage, A ~ 1/t. In the long-time limit, 
A = (40/9)(fcbox/fcflow) 2 <C 1. The spectral exponent at 
k S> 1 is given by Eq. (JJJ and is only slightly shallower 
than — 1. The spectrum at small scales is, thus, very 
similar to the Batchelor spectrum for the forced scalar 
turbulence H3- The physical reason for this similarity is 
that the modes at the flow scale and above act as a slowly 
decaying source to the small-scale part of the spectrum, 
thus making the small-scale physics similar to the forced 
case (cf. Q). 

C. Decay of scalar variance (C2 > 0) 

Finally, we estimate the decay of the total scalar vari- 
ance. The total variance is the integral of the wave- 
number spectrum, which is made up of the nonsingu- 
lar spectrum [Eq. H21|) ] and the two singular peaks. To 
obtain the contribution of the latter to the total vari- 
ance, we must take into account their time-dependent 
widths. We see that the contributions from the flow- 
scale band and from the nonsingular spectrum are al- 
ways of the same order ~ A(t)Ti(t) [taking into account 
that the spectrum at k > 1 is f(k) ~ fc- 4 +( 3 / 4 ) A U p to 
k n ~ Pe 1 / 2 gives an extra factor of Pe^ 3 / 8 ^ to the non- 
singular contribution]. Therefore, during the transient 
stage (t <C l/^box)' tms P ar t °f the variance decays as 
1/t 5 / 2 . In the long-time limit (t > l/fc^ ox ), we have 
AO)Ti(i) - fcboxTiO) oc t- 1 exp(-At). The large-scale 
peak always decays as tT\(t). In the transient stage, 
its width is ~ 1/vr, which means that its contribution 
~ 1 /t 3 / 2 dominates the rest of the spectrum by a factor 
of ~ t and determines the overall decay law of the total 
variance [35l |. In the long-time limit, the width of the box 
mode is ~ l/khoxt, so its decay law is ~ Ti(*)/fcbcoc> i- e -j 
it has the same time dependence as the rest of the spec- 
trum, but its contribution to the total variance exceeds 
that of all other modes by a factor of ~ 1/fcj^. 

Note that these arguments can be tested for consis- 
tency with the conservation law for the scalar variance in 
the following way. Pick some wave number 1 <C k <C k n . 
Because of the dominance of the box mode, the scalar 
variance integrated up to k is the same as the total vari- 
ance. Its time derivative must be equal to the flux of vari- 
ance through k [the negative of the expression in square 
brackets in Eq. (J1J] because dissipation is negligible at 
k <^k n . It is easy to check that this, indeed, holds true. 
This argument emphasizes that the turbulent diffusion at 
the box scale, which we have shown to control scalar de- 
cay, is not a dissipation mechanism, but rather describes 
transfer of scalar variance to small scales, where all the 
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FIG. 2: Evolution of normalized spectrum T(t,k)/E(t) 
for fcbox = 0.1. The initial spectrum was T(0, k) = 
k 4 exp[— (k/k v ) 2 ]. The long-term solution for fcbox = 0.01 
is also shown. Spectra at k < 1 are steeper than k 1 due to 
higher-order corrections: they can be well fitted by polyno- 
mials in the form k + ak 2 + . . . [Eq. The blip at k = 3 
is the point where the solutions of Eq. JHJ at k < 3 and of 
Eq. Q at k > 3 are spliced together. This device allows us to 
achieve higher Pe. We have checked that moving around the 
splicing wave number or solving the full integral equation in 
the entire domain (at lower Pe) do not change the solution. 



dissipation is done by molecular diffusion. 



D. Case of C7 2 = 

When the Corrsin invariant is zero, the theory devel- 
oped above has to be modified somewhat. The dominant 
term in the small-fc solution is ~ fc 4 and includes con- 
tributions both from the initial spectrum and from the 
initial finite amount of nonlinear transfer of scalar vari- 
ance from the flow-scale mode: Equation l|12f> is replaced 
by 

T(t,k) = [C* 4 + $(i,fc)]fc 4 e- fc2 *, (23) 
$(t,fc) = ~ [ dt'e k2t ' f dz(l~ z 2 )T{t',l + kz). (24) 

In the transient stage (1 « f C V^box)' we nave i at 
k < $(t,k) ~ j*dt'Ti(t') -> const as long as 

Ti(t) decays faster than 1/t [cf. Eq. ifTijl], Since 
also remains finite for fc ~ 1/ y/i, we can estimate the flow 
scale modes T(t, 1 + q) by substituting the solution lll'-'il) 
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into Eq. (|llfl and assuming that <I>(i, fc) is approximately 
constant at the values of A; that matter. This will produce 
errors of order unity in the prefactors but yield correct 
scalings and asymptotic decay laws. Note that because 
the nonlinear-transfer contribution is ~ fc 4 , the higher- 
order terms in the Taylor expansion for the initial spec- 
trum do not affect the asymptotic solutions. 

Further derivation for the Ci = case follows the same 
general scheme as the theory for C2 > 0. In the transient 
stage, Eq. ijT3|) is replaced by 

T{t,l + q)~ 2 -tjfie-f\ (25) 

where we have dropped prefactors of order unity. Note 
that the transient-stage decay of the flow-scale peak, 
Ti(t) ~ 1/t 3 , is faster then for the C 2 > 0. This is 
the only important change that results from the vanish- 
ing of C2. The intermediate asymptotic (1/y/t <C k <C 1) 
is again T(t,k) ~ T 1 (t)t- 1 / 2 k [cf. Eq. fT5fl. In the long- 
time limit, Eq. (|l(j[l is replaced by 

T(t,l + q)-'^(kl ox - q 2 )e- k ^, (26) 
the box mode decays according to [cf. Eq. (|17|l ] 

T(t, fc box ) ~ fcg„e-*S«« [C 4 + O(k£ ox In t)] , (27) 

and the intermediate asymptotic at fcb ox "C k <C 1 is 
T(t, fc) ~ Ti(t)fc box fc [cf. Eq. ©]. 

The developments for the nonsingular spectrum and 
for the total scalar variance are exactly as described in 
Sec. ITTTBl and Sec. ImCl The only change is in the 
transient-stage decay laws due to faster decay of the flow- 
scale mode [Eq. <|25[) ]: the contribution to the scalar vari- 
ance from the flow scale and the nonsingular spectrum 
is ~ 1/t 7 / 2 ; the contribution from the small-fc peak is 
~ 1/t 5 / 2 . The latter dominates and, therefore, deter- 
mines the decay law for the total variance. The behavior 
in the long-time limit is the same as for C2 > 0. 

E. Numerical solution 

We have checked our analytical solution by solving 
Eq. (JSJ) numerically. The powerlike decay laws for the 
flow-scale mode T(t, 1) and for the total variance are in 
good agreement with theory, but they can only can be 
seen if the scale separation between the box and the flow 
is sufficiently large [fcbox ~ 0.01; see Fig. ^a)]. Spec- 
tra (Fig. [5J and long-term decay rates [Fig. db)] are 
clearly in agreement with theory already at moderate 
scale separations. A simple way to estimate the thresh- 
old at which the decay rate ceases to be determined by 
box-scale diffusion is by requiring that the box decay rate 
should be smaller than the Lagrangian-stretching value: 
A = (40/9)fc^ ox < 1, so fcbox < 0.47. Consistent with 
this estimate, the box value still works for fcbox = 1/3 
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FIG. 3: Structure of the "strange mode." The spectrum is 
from a run with fcbox = 0.1, k v = 100 (the same as in Fig. [5] 
t — 1000). Gray arrows show directions of couplings (transfer 
of scalar variance) . 

[Fig- DIb); cf. 17]]. Obviously, Eq. © itself with the 
cutoff at fcbox is only technically valid when fcbox <C 1, 
but we see that it continues to yield reasonable solutions 
even at moderate fcbox- 



IV. DISCUSSION 

We now have a physical picture of the "strange mode" : 
the small-fc peak (singularity at the box scale) serves as a 
slowly decaying source to the flow-scale mode (singularity 
at fc = 1), which, in turn, is mixed by the random flow 
and thus excites the nonsingular modes at small [Eq. JSJ] 
and large [Eq. (|22|l ] scales. The structure of the spectrum 
is illustrated in Fig. |3J 

The low-wave-number behavior of the decaying scalar 
field was previously analyzed in a heuristic way by Ker- 
stein and McMurtry [24| (see also [25| for a treatment 
based on one of the turbulence closure schemes, which 
gives mostly similar results). They considered advection 
by a narrow-band (i.e., single-scale) forced random flow 
in an unbounded domain — i.e., in the regime that we call 
the transient (powerlike-decay) stage. They recognized 
the defining role of coupling between the large scales 
(k < fcpoak) and the flow scale (fc = fcfl ow = 1) and derived 
the fc 4 scaling at fc <g; fc pca k with a exponential fall off at 
fc > fcpcak [Eq. (|12fl ] and the ensuing powerlike-decay laws 
for the case of C2 = (see the end of Sec. IIII Dfl . For 
the intermediate rage fc p0 ak *C k <C 1, they predicted 
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a fc 2 spectrum (in 3D) — in contrast to our fc 1 result 
[Eqs. G5J and @ and Sec. HTTP] . The reason for the 
discrepancy is as follows. The analysis of |24j is based on 
Taylor-expanding the flow around k — 1 — i.e., in terms 
of our theory — setting T(t, 1 + q) ~ T(t, 1) in Eq. J5J, 
which gives S(t,k) ~ k 4 T(t, 1). If we had used the re- 
sulting equation to solve for T(f, fc) at fc pca k <C fc <C 1, we 
would also have obtained T(t, fc) ~ A: 2 . However, as we 
have seen above, the width of the flow-scale singularity is 
~ frpeak [Eqs. (|T3|> and so Taylor expansion cannot 

be used in Eq. © for k ^> fc pcak . In this intermediate 
range, the integral in Eq. 10 must instead be replaced 
by the integral over the entire flow-scale peak, resulting 
in our fc 1 scaling. The k 2 term enters as a correction due 
to the interaction between nonsingular modes [Eq. Ij22(l ]. 

Finally, let us comment on our modeling assumptions. 
The white- noise approximation might appear drastic: the 
correlation time of any realistic flow is comparable to 
the flow time scale ~ (ufcfl ow ) _1 ~ n^ 1 - However, since 
the scalar decay time is much longer than the flow time 
scale (provided fcb ox <C fcfl w)j the white-noise model ap- 
pears reasonable. We believe it also correctly captures 
the small-scale structure: the key factor here is the statis- 
tics of fluid displacements, which are integrals of velocity 
and are finite-time correlated even for a white-in-time 
velocity. 

Our model flow was single scale. Although such flows 
can be set up in the laboratory 0, (3(| , the real- world 
mixing problems usually contain (at sufficiently small 
scales) a wide (inertial) scale range of three-dimensional 



turbulent motions. While the variance spectrum in the 
inertial range should follow the Obukhov-Corrsin fc~ 5 / 3 
law jlEISSl and there will be another transient power like- 
decay stage ^3113; HHHIiHiii the long-term decay (after 
the scalar variance reaches k < fcfl ow ) should still be qual- 
itatively described by our theory. Another modification 
that results from the relaxation of the single-scale as- 
sumption concerns the intermediate wave-number range 
fcpcak < < fcflow As noted in [24| and confirmed in 
pipe-flow mixing experiments |3lT |. the interaction be- 
tween the fc = fcpcak mode and the low-wave-number tail 
of the kinetic- energy spectrum (~ fc 4 ) can change the 
scaling of the scalar-variance spectrum in this range. 

In conclusion, we emphasize that, in any laboratory 
experiment aiming to test our results, the stirring must 
be done at scales substantially smaller than the system 
size to ensure that fcbox "C fcflow It was just such a set up 
(in 2D) that allowed Voth et al. [13 to show experimen- 
tally that the global mixing rate was much smaller than 
that predicted by the Lagrangian-stretching theories and 
consistent with the box-scale turbulent-diffusion rate — 
precisely the point the theory presented above is meant 
to demonstrate. 
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